Comparison of Ethos template-based planning and AI-based dose prediction: General performance, patient optimality, and limitations.

PURPOSE: Ethos proposes a template-based automatic dose planning (Etb) for online adaptive radiotherapy. This study evaluates the general performance of Etb for prostate cancer, as well as the ability to generate patient-optimal plans, by comparing it with another state-of-the-art automatic planning method, i.e., deep learning dose prediction followed by dose mimicking (DP + DM). MATERIALS: General performances and capability to produce patient-optimal plan were investigated through two studies: Study-S1 generated plans for 45 patients using our initial Ethos clinical goals template (EG_init), and compared them to manually generated plans (MG). For study-S2, 10 patients which showed poor performances at study-S1 were selected. S2 compared the quality of plans generated with four different methods: 1) Ethos initial template (EG_init_selected), 2) Ethos updated template-based on S1 results (EG_upd_selected), 3) DP + DM, and 4) MG plans. RESULTS: EG_init plans showed satisfactory performance for dose level above 50 Gy: reported mean metrics differences (EG_init minus MG) never exceeded 0.6 %. However, lower dose levels showed loosely optimized metrics, mean differences for V30Gy to rectum and V20Gy to anal canal were of 6.6 % and 13.0 %. EG_init_selected showed amplified differences in V30Gy to rectum and V20Gy to anal canal: 8.5 % and 16.9 %, respectively. These dropped to 5.7 % and 11.5 % for EG_upd_selected plans but strongly increased V60Gy to rectum for 2 patients. DP + DM plans achieved differences of 3.4 % and 4.6 % without compromising any V60Gy. CONCLUSION: General performances of Etb were satisfactory. However, optimizing with template of goals might be limiting for some complex cases. Over our test patients, DP + DM outperformed the Etb approach.

Dosimetrically triggered adaptive radiotherapy for head and neck cancer: Considerations for the implementation of clinical protocols.

PURPOSE: Defining dosimetric rules to automatically detect patients requiring adaptive radiotherapy (ART) is not straightforward, and most centres perform ad-hoc ART with no specific protocol. This study aims to propose and analyse different steps to design a protocol for dosimetrically triggered ART of head and neck (H&N) cancer. As a proof-of-concept, the designed protocol was applied to patients treated in TomoTherapy units, using their available software for daily MVCT image and dose accumulation. METHODS: An initial protocol was designed by a multidisciplinary team, with a set of flagging criteria based only on dose-volume metrics, including two action levels: (1) surveillance (orange flag), and (2) immediate verification (red flag). This protocol was adapted to the clinical needs following an iterative process. First, the protocol was applied to 38 H&N patients with daily imaging. Automatic software generated the daily contours, recomputed the daily dose and flagged the dosimetric differences with respect to the planning dose. Second, these results were compared, by a sensitivity/specificity test, to the answers of a physician. Third, the physician, supported by the multidisciplinary team, performed a self-analysis of the provided answers and translated them into mathematical rules in order to upgrade the protocol. The upgraded protocol was applied to different definitions of the target volume (i.e. deformed CTV + 0, 2 and 4 mm), in order to quantify how the number of flags decreases when reducing the CTV-to-PTV margin. RESULTS: The sensitivity of the initial protocol was very low, specifically for the orange flags. The best values were 0.84 for red and 0.15 for orange flags. After the review and upgrade process, the sensitivity of the upgraded protocol increased to 0.96 for red and 0.84 for orange flags. The number of patients flagged per week with the final (upgraded) protocol decreased in median by 26% and 18% for red and orange flags, respectively, when reducing the CTV-to-PTV margin from 4 to 2 mm. This resulted in only one patient flagged at the last fraction for both red and orange flags. CONCLUSION: Our results demonstrate the value of iterative protocol design with retrospective data, and shows the feasibility of automatically-triggered ART using simple dosimetric rules to mimic the physician's decisions. Using a proper target volume definition is important and influences the flagging rate, particularly when decreasing the CTV-to-PTV margin.

Automated clinical decision support system with deep learning dose prediction and NTCP models to evaluate treatment complications in patients with esophageal cancer.

BACKGROUND AND PURPOSE: This study aims to investigate how accurate our deep learning (DL) dose prediction models for intensity modulated radiotherapy (IMRT) and pencil beam scanning (PBS) treatments, when chained with normal tissue complication probability (NTCP) models, are at identifying esophageal cancer patients who are at high risk of toxicity and should be switched to proton therapy (PT). MATERIALS AND METHODS: Two U-Net were created, for photon (XT) and proton (PT) plans, respectively. To estimate the dose distribution for each patient, they were trained on a database of 40 uniformly planned patients using cross validation and a circulating test set. These models were combined with a NTCP model for postoperative pulmonary complications. The NTCP model used the mean lung dose, age, histology type, and body mass index as predicting variables. The treatment choice is then done by using a ΔNTCP threshold between XT and PT plans. Patients with ΔNTCP ≥ 10% were referred to PT. RESULTS: Our DL models succeed in predicting dose distributions with a mean error on the mean dose to the lungs (MLD) of 1.14 ± 0.93% for XT and 0.66 ± 0.48% for PT. The complete automated workflow (DL chained with NTCP) achieved 100% accuracy in patient referral. The average residual (ΔNTCP ground truth - ΔNTCP predicted) is 1.43 ± 1.49%. CONCLUSION: This study evaluates our DL dose prediction models in a broader patient referral context and demonstrates their ability to support clinical decisions.

Linear Energy Transfer Incorporated Spot-Scanning Proton Arc Therapy Optimization: A Feasibility Study.

PURPOSE: To integrate dose-averaged linear energy transfer (LETd) into spot-scanning proton arc therapy (SPArc) optimization and to explore its feasibility and potential clinical benefits. METHODS: An open-source proton planning platform (OpenREGGUI) has been modified to incorporate LETd into optimization for both SPArc and multi-beam intensity-modulated proton therapy (IMPT) treatment planning. SPArc and multi-beam IMPT plans with different beam configurations for a prostate patient were generated to investigate the feasibility of LETd-based optimization using SPArc in terms of spatial LETd distribution and plan delivery efficiency. One liver and one brain case were studied to further evaluate the advantages of SPArc over multi-beam IMPT. RESULTS: With similar dose distributions, the efficacy of spatially optimizing LETd distributions improves with increasing number of beams. Compared with multi-beam IMPT plans, SPArc plans show substantial improvement in LETd distributions while maintaining similar delivery efficiency. Specifically, for the liver case, the average LETd in the GTV was increased by 124% for the SPArc plan, and only 9.6% for the 2-beam IMPT plan compared with the 2-beam non-LETd optimized IMPT plan. In case of LET optimization for the brain case, the SPArc plan could effectively increase the average LETd in the CTV and decrease the values in the critical structures while smaller improvement was observed in 3-beam IMPT plans. CONCLUSION: This work demonstrates the feasibility and significant advantages of using SPArc for LETd-based optimization, which could maximize the LETd distribution wherever is desired inside the target and averts the high LETd away from the adjacent critical organs-at-risk.

Deep learning dose prediction for IMRT of esophageal cancer: The effect of data quality and quantity on model performance.

PURPOSE: To investigate the effect of data quality and quantity on the performance of deep learning (DL) models, for dose prediction of intensity-modulated radiotherapy (IMRT) of esophageal cancer. MATERIAL AND METHODS: Two databases were used: a variable database (VarDB) with 56 clinical cases extracted retrospectively, including user-dependent variability in delineation and planning, different machines and beam configurations; and a homogenized database (HomDB), created to reduce this variability by re-contouring and re-planning all patients with a fixed class-solution protocol. Experiment 1 analysed the user-dependent variability, using 26 patients planned with the same machine and beam setup (E26-VarDB versus E26-HomDB). Experiment 2 increased the training set by groups of 10 patients (E16, E26, E36, E46, and E56) for both databases. Model evaluation metrics were the mean absolute error (MAE) for selected dose-volume metrics and the global MAE for all body voxels. RESULTS: For Experiment 1, E26-HomDB reduced the MAE for the considered dose-volume metrics compared to E26-VarDB (e.g. reduction of 0.2 Gy for D95-PTV, 1.2 Gy for Dmean-heart or 3.3% for V5-lungs). For Experiment 2, increasing the database size slightly improved performance for HomDB models (e.g. decrease in global MAE of 0.13 Gy for E56-HomDB versus E26-HomDB), but increased the error for the VarDB models (e.g. increase in global MAE of 0.20 Gy for E56-VarDB versus E26-VarDB). CONCLUSION: A small database may suffice to obtain good DL prediction performance, provided that homogenous training data is used. Data variability reduces the performance of DL models, which is further pronounced when increasing the training set.

Dose prediction with deep learning for prostate cancer radiation therapy: Model adaptation to different treatment planning practices.

PURPOSE: This work aims to study the generalizability of a pre-developed deep learning (DL) dose prediction model for volumetric modulated arc therapy (VMAT) for prostate cancer and to adapt the model, via transfer learning with minimal input data, to three different internal treatment planning styles and one external institution planning style. METHODS: We built the source model with planning data from 108 patients previously treated with VMAT for prostate cancer. For the transfer learning, we selected patient cases planned with three different styles, 14-29 cases per style, in the same institution and 20 cases treated in a different institution to adapt the source model to four target models in total. We compared the dose distributions predicted by the source model and the target models with the corresponding clinical plan dose used for patient treatments and quantified the improvement in the prediction quality for the target models over the source model using the Dice similarity coefficients (DSC) of 0% to 100% isodose volumes and the dose-volume-histogram (DVH) parameters of the planning target volume and the organs-at-risk. RESULTS: The source model accurately predicts dose distributions for plans generated in the same source style, but performs sub-optimally for the three different internal and one external target styles, with the mean DSC ranging between 0.81-0.94 and 0.82-0.91 for the internal and the external styles, respectively. With transfer learning, the target model predictions improved the mean DSC to 0.88-0.95 and 0.92-0.96 for the internal and the external styles, respectively. Target model predictions significantly improved the accuracy of the DVH parameter predictions to within 1.6%. CONCLUSION: We demonstrated the problem of model generalizability for DL-based dose prediction and the feasibility of using transfer learning to solve this problem. With 14-29 cases per style, we successfully adapted the source model into several different practice styles. This indicates a realistic way forward to widespread clinical implementation of DL-based dose prediction.

Dosimetric evaluation of synthetic CT generated with GANs for MRI-only proton therapy treatment planning of brain tumors.

PURPOSE: The purpose of this study was to address the dosimetric accuracy of synthetic computed tomography (sCT) images of patients with brain tumor generated using a modified generative adversarial network (GAN) method, for their use in magnetic resonance imaging (MRI)-only treatment planning for proton therapy. METHODS: Dose volume histogram (DVH) analysis was performed on CT and sCT images of patients with brain tumor for plans generated for intensity-modulated proton therapy (IMPT). All plans were robustly optimized using a commercially available treatment planning system (RayStation, from RaySearch Laboratories) and standard robust parameters reported in the literature. The IMPT plan was then used to compute the dose on CT and sCT images for dosimetric comparison, using RayStation analytical (pencil beam) dose algorithm. We used a second, independent Monte Carlo dose calculation engine to recompute the dose on both CT and sCT images to ensure a proper analysis of the dosimetric accuracy of the sCT images. RESULTS: The results extracted from RayStation showed excellent agreement for most DVH metrics computed on the CT and sCT for the nominal case, with a mean absolute difference below 0.5% (0.3 Gy) of the prescription dose for the clinical target volume (CTV) and below 2% (1.2 Gy) for the organs at risk (OARs) considered. This demonstrates a high dosimetric accuracy for the generated sCT images, especially in the target volume. The metrics obtained from the Monte Carlo doses mostly agreed with the values extracted from RayStation for the nominal and worst-case scenarios (mean difference below 3%). CONCLUSIONS: This work demonstrated the feasibility of using sCT generated with a GAN-based deep learning method for MRI-only treatment planning of patients with brain tumor in intensity-modulated proton therapy.

Accelerated robust optimization algorithm for proton therapy treatment planning.

PURPOSE: Robust optimization is a computational expensive process resulting in long plan computation times. This issue is especially critical for moving targets as these cases need a large number of uncertainty scenarios to robustly optimize their treatment plans. In this study, we propose a novel worst-case robust optimization algorithm, called dynamic minimax, that accelerates the conventional minimax optimization. Dynamic minimax optimization aims at speeding up the plan optimization process by decreasing the number of evaluated scenarios in the optimization. METHODS: For a given pool of scenarios (e.g., 63 = 7 setup  × 3 range  × 3 breathing phases), the proposed dynamic minimax algorithm only considers a reduced number of candidate-worst scenarios, selected from the full 63 scenario set. These scenarios are updated throughout the optimization by randomly sampling new scenarios according to a hidden variable P, called the "probability acceptance function," which associates with each scenario the probability of it being selected as the worst case. By doing so, the algorithm favors scenarios that are mostly "active," that is, frequently evaluated as the worst case. Additionally, unconsidered scenarios have the possibility to be re-considered, later on in the optimization, depending on the convergence towards a particular solution. The proposed algorithm was implemented in the open-source robust optimizer MIROpt and tested for six four-dimensional (4D) IMPT lung tumor patients with various tumor sizes and motions. Treatment plans were evaluated by performing comprehensive robustness tests (simulating range errors, systematic setup errors, and breathing motion) using the open-source Monte Carlo dose engine MCsquare. RESULTS: The dynamic minimax algorithm achieved an optimization time gain of 84%, on average. The dynamic minimax optimization results in a significantly noisier optimization process due to the fact that more scenarios are accessed in the optimization. However, the increased noise level does not harm the final quality of the plan. In fact, the plan quality is similar between dynamic and conventional minimax optimization with regard to target coverage and normal tissue sparing: on average, the difference in worst-case D95 is 0.2 Gy and the difference in mean lung dose and mean heart dose is 0.4 and 0.1 Gy, respectively (evaluated in the nominal scenario). CONCLUSIONS: The proposed worst-case 4D robust optimization algorithm achieves a significant optimization time gain of 84%, without compromising target coverage or normal tissue sparing.